Spin gap of the three-leg S* = 3/2 Heisenberg tube 
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The ground-state properties of the three- leg 5* = 3/2 Heisenberg tube are studied using the 
density-matrix renormalization group method. We find that the spin-excitation gap associated with 
a spontaneous dimerization opens for the whole coupling regime, as seen in the three-leg 5* = 1/2 
Heisenberg tube. However, in contrast to the case ol S — \ tube, the gap increases very slowly 
with increasing the rung coupling and its size is only a few % or less of the leg exchange interaction 
in the weak- and intermediate-coupling regimes. We thus argue that, unless the rung coupling is 
substantially larger than the leg coupling, the gap may be quite hard to be observed experimentally. 
We also calcuate the quantized Berry phase to show that there exist three kinds of valence-bond-solid 
states depending on the ratio of leg and rung couplings. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Cx, 75.40.Mg 



For a long time, the exotic phenomena emerged 
from geometrical frustration have been fascinating but 
challenging subjects of research in condensed matter 
physicsii The peculiar dilemma of frustrated systems 
generally comes from a highly-degenerate ground state 
in the classical sense. Here, we know that to resolve 
it comes essentially back to how the degeneracy is re- 
moved or how the frustration is minimized by taking the 
quantum fluctuations into account. The simplest exam- 
ple for the spin frustration may be the 120° structure of 
antiferromagnetic triangle. In the context of triangular- 
lattice 5 = 1/2 antiferromagnet a spin-liquid state was 
proposed by Andcrsoni^ As a related issue, odd-leg spin 
tube systems such as Na2V307 (Ref. d) and [(CuCl2 
tachH)3Cl]Cl2 (Ref. 3) have attracted much attention for 
the last few years. One could say that odd- leg spin lad- 
ders belong to the same universality class as single chain 
does; thus, the ground state is comprehended as a gap- 
less spin liquid or a Tomonaga-Luttinger (TL) liquid*^ 
However, if the periodic boundary conditions (PBC) are 
applied in the rung direction, i.e., a tube is shaped, the 
spin states are dramatically changed due to spin frustra- 
tion in the polygonal ring with odd number of rungs. 

Quite recently, the hexagonal compound CsCrF4 
(Ref. 0), which is an ideal three- leg spin tube system 
formed by Cr'^^ ions, has been reexamined experimen- 
tally from the point of view of spin frustrationi^ Since 
the magnetic moment comes from the state of Cr'^^ 
ion, the magnitude of spin on each site is S" = | . By per- 
forming magnetic susceptibility, heat capacity C(T), and 
electron spin resonance measurements, it was reported 
that no magnetic long-range order is observed down to 
T = 1.3K. In particular, a gapless spin-liquid state (or 
a TL liquid state) is indicated from the finite T-linear 
component of C{T); this result raises a more absorbing 
question because a gapped ground state is expected in 
odd-leg spin-half-integer spin tube system<^ The need for 
an investigation of odd-leg spin tube system with S = § 




FIG. 1. Lattice structure of three-leg Heisenberg tube. 



is therefore obvious in order to figure out this puzzle. 

In this paper, we thus consider three-leg 5" = | Heisen- 
berg tube. The Hamiltonian is given as 



H — J\\ Sa,i ■ Sa,i+1 + J± Sa.i ■ Sa'.i, 

a=l i i a(^a') 

where Sa,i is a spin-| operator at leg a(= 1,2,3) and 
rung i. J|| and Jj_ are the nearest-neighbor exchange 
interactions in the leg and rung directions, respectively 
(see FigH]). We take J\\ = 1 as the unit of energy here- 
after. In order to investigate the ground-state proper- 
ties of the system ([T]) , the density- matrix renormalization 
group (DMRG) techniquo^i is employed. As necessary, 
the PBC or the open-end boundary conditions (OBC) 
are chosen in the leg direction. Using the OBC (PBC), 
we study the tubes with several kinds of length L = 12 
to 48 (i = 8 to 24) keeping m = 1200 to 2600 {m = 1600 
to 3200) density-matrix eigenstates in the renormaliza- 
tion procedure; in this way, the typical truncation er- 
ror, i.e., the discarded weight, is 2 x 10^*^ — 1 x 10~'^ 
(3 X 10"^ - 7 X 10~^). We note that the system length 
must be even and is better to be kept in L = 4/ or 4Z + 2 
with Z=intcger for systematic extrapolation of calculated 
quantities into the thermodynamic limit. Moreover, an 
extrapolation to m — > cxd for each calculation is necessary 
because the DMRG trial wave function slowly approaches 
the exact one with increasing m due to the large degrees 
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FIG. 2. Schematic pictures of the VBS configurations. Each 
small solid circle and connection with line denote a spin-^ 
variable and a singlet pair, respectively. The large open circle 
represents a spin-| operation which symmetrizes three spin-^ 
variables. The configurations (I) and (II) are VBS states in 
small and large bond-alternation regimes of the ID S = | 
Heisenberg model, respectively. The configurations (III)-(V) 
are possible candidates of VBS states in the three-leg 5=1 
Heisenberg tube (details are given in the text and Fig. [5]). 

of freedom ^ 4^^ and strong spin frustration in our sys- 
tem. All the calculated quantities in this paper are ex- 
trapolated to the limit m — > oo; thus, for example, the 
maximum error in the ground-state energy is estimated 
to be less than 1 x lO^'^. 

The first thing wc think of when considering the sys- 
tem ([T]) might be the topological similarity to the other 
odd-leg half-integer-spin Heisenberg tubes. The sim- 
plest case, i.e., three- leg S = ^ Heisenberg tube, has 
been well-studied, and the ground state is known to be 
gapped where the system is spontaneously dimerized in 
the leg direction to relax the intra-rung spin frustra- 
tion.— liSrJ^ This can be naturally understood by analogy 
with the gap-opening mechanism in the one-dimensional 
(ID) 5* = i spin-Peierls Heisenberg modelj^^ Hence, in 
the case of 5 = | as well, it would be best to start with 
bond-alternated single chain problem, namely, the ID 
5 = 1 spin-Peierls Heisenberg model. The Hamiltonian 

is written as H — ^f^i[l — {~iyS]Si ■ S'i+i where Si 
is a spin-| operator at site i and i5(> 0) is the strength 
of bond alternation. The low-energy physics of this sys- 
tem has been fundamentally elucidated J^ii^ Across the 
critical point 5 « 0.42, two kinds of valence-bond-solid 
(VBS) phases appear in the ground state; for the larger 
alternation (5 > 0.42), the VBS state is essentially writ- 
ten as a direct product of 'simple' spin-Peierls singlet 
bonds [Fig. HJH)], whereas for the smaller alternation 
{S < 0.42), it is expressed as a combined state of the 
spin-Peierls singlet and S" = 1 Haldane-like-gapped con- 
figurations [Fig. [2Il)]. And, the ground state is always 
gapped except at the critical point. Now therefore, get- 
ting back to our system ([T} , if the spontaneous dimeriza- 
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FIG. 3. Dimerization order parameter 7? as a function of 
J±. Inset: Semilog plot of D vs. 1/J±. D saturates to 0.2764 
in the large Jx limit. The data for small Jx is fitted by a 
function D = aexp(— /3/Jx) with a = 0.25 and jS = 13.2. 

tion occurs as in the three-leg S = \ Heisenberg tube, it 
is likely that a gapped ground state is also obtained here. 

Then, we will simply evaluate a dimerization order pa- 
rameter to check the presence or absence of long-ranged 
spin-Peierls ordering in our system. Since a spin-Peierls 
state is characterized as an alternating formation of spin- 
singlet pairs in the leg direction, we focus on the nearest- 
neighbor spin-spin correlations, 

5(*) = -(5'a,,-5'„,,+i), (2) 

where (■ • ■) denotes an expectation value in the ground 
state. Note that this quantity is independent of a. With 
applying the OBC, the translational symmetry is broken 
due to the Friedel oscillation and the spin-Peierls state 
is directly observable as a ground state. In general, the 
Friedel oscillation in the center of the system decays as 
a function of the system length. If the amplitude at the 
center of the system 

D{L)^\S{L/2)-S{L/2 + l)\ (3) 

persists for arbitrarily long system length, it corresponds 
to a long-ranged dimerization order which indicates the 
spin-Pcicrls ground state. The order parameter is thus 
defined as an extrapolated value into the thermodynamic 
limit, 

D = lim D{L). (4) 

L— s-oo 

In Fig. [3l the extrapolated values D are plotted as 
a function of Jj_. We see that it increases gradually 
at J_L ^ 5, almost linearly at 5 ^ J_l ^ 15, and then 
go into saturation at Z? = 0.2764 in the large J± limit. 
These different behaviors could be interpreted in terms 
of different VBS state for each the J± regime, as in the 
S = ^ Heisenberg tubei^i This point will be clarified be- 
low by examining the Berry phase. A remaining question 
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FIG. 4. (a) Finite-size scaling of A(I/) as a function of l/L. 
The lines are the polynomial fittings, (b) Extrapolated val- 
ues of A(I/) to the thermodynamic limit 1/L — > 0. Inset: 
extended figure for < Jj. < 4. The error bars give differ- 
ences between the second-order and cubic polynomial fittings 
for A(L). 



would be whether the order parameter remains finite for 
small J_L regime (we have not successfully obtained D for 
Jj^ < 3 due to its smallness). We find that the behavior 
seems just like the Berezinskii-Kosterlitz-Thouless type 
transition, D = 0.25 exp(— 13.2/ Jj^), from the single log- 
arithmic plot (see the inset of Fig. ; it may imply that 
the order parameter is exponentially small but finite at 
< J_i < 3. We conclude that in the wide range of J_l 
the dimerization order occurs and the spin-excitation gap 
is expected to be finite there. 

Let us then estimate the spin-excitation gap. Of par- 
ticular interest is the evolution of the gap onto the ratio 
between leg and rung couplings. The gap is defined as 



A = £;i(L)-So(i), 



A= lim A(L), 

L— ^CX3 



(5) 



where Eo{L) and Ei{L) are energies of the ground state 
{S = 0) and first triplet excited state [S = 1) for the 
system with length L, respectively. In Fig. HJa), we 
plot the system-size dependence of the spin-excitation 
gap calculated with the OBC in full circles. We see 
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FIG. 5. Schematic phase diagram of the three-leg S = ^ 
Heisenberg tube, classified by the Berry phases on the leg 
bond (7icg) and rung bond (7iGg). The Roman numbers cor- 
respond to the VBS states shown in Fig. [21 



that the values of A(L) can be smoothly extrapolated 
to the thermodynamic limit 1/L oo. The extrapo- 
lated values A. using a cubic polynomial extrapolation 
for A(L), are shown in Fig. SKb) as a function of J±. 
As expected, J_L-dependence of A looks rather similar to 
that of D] namely, it increases slowly at 5, rapidly 

at 5 ^ ^ 15, and then saturates to A = 0.6661 in 
the strong-coupling limit J_i = cxd. This is because the 
spin-excitation gap is essentially equivalent to a binding 
energy of most weakly bounded spin-singlet pair in the 
system and it is approximately scaled with the dimeriza- 
tion strength. 

Here, it is very instructive to compare the gap with 
that of the S ~ ^ tube, which is also shown in Fig. Hl^b). 
Two remarks are made by the comparison: (i) Although 
it may be rather natural, the gap in the strong-coupling 
limit seems to be scaled with the magnitude of spin, 
A(J^ oo) cx S, (ii) the gap for = | increases much 
more slowly with increasing Jj^ in the weak-coupling 
regime. As a result, only a few % of the leg exchange 
interaction remains even at J±/ J|| = 5. Hence, we argue 
that, unless the ratio J_l/./|| is sizably large, it may be 
difficult to detect the gap experimentally. For CsCrF4, 
the leg exchange interaction is estimated to be a few 10 
to 100 K by comparing the experimental peak position of 
magnetic susceptibility to numerical analysis for the ID 

= I Heisenberg chaiufii and the gap is only a few K 
at the outside even for J±/ J\\ =5. 

As described above, we obtain the gapped ground state 
for the whole J± region with applying the OBC. It would 
mean that our system never contain the Haldane-type 
VBS state [Fig. E] (I)] in the three chains. This is be- 
cause the gap cannot open due to free edge spins created 
with the OBC if the state (I) is included. But to be 
sure, we shall confirm it by estimating the gap with the 
PBC. The obtained results are shown with open squares 
in FigHl We see that the extrapolated values are in good 
agreement to those with the OBC and it is confirmed 
that the state (I) does not exist at any VBS state in our 
system. Now, it is a fair question then to ask which kind 
of VBS state is formed and how it changes with varying 
J±. 

Finally, we investigate the quantized Berry phase for 
determining topological configuration of VBS ground 
state. The Berry phase is defined by 



-I I A{(P)d<p, 





(6) 
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where is the Abehan Berry connection, A{cj)) = 

{tp^\d^ip^) with the ground state The Berry phase 

is quantized as or tt (mod 27r) if the system has spin 
gap during the adiabatic continuation and time rever- 
sal symmetry; and "undefined" if a gapless excitation 
exists. We introduce a local perturbation by a local 
twist of the nearest-neighbor connection, Sa i ■ Sa' j 
U^-''^Sl,S-,^^+e^'^S-X',j) + Sl.S^^,,y In this paper, 
a couple of clusters with L = 2 and 4 are analyzed for 
this quantity. A dimerized pair of triangles, including six 
spins, from the clusters arc picked up (it is the cluster 
itself for L = 2), and the Berry phases of the leg bond 
(7ieg) for a = a' , j = i + 1 and of the rung bond (7rung) 
for a ^ a',j = i are evaluated. We call the spin-singlet 
pair on the leg (rung) bond "on-leg (on-rung) pair" . 

The Berry phases are obtained as follows: 
(7iog,7rung) = (7r,0) at < J± < 1 (0 < J± < 0.5), 
(7iog,7rung) = (7r,7r) at 1 < J_L < 15.3 (1 < J_l < 18), 
and (7iog,7rung) = (0,7r) at > 15.3 {J±_ > 18) 
for i = 2 (L = 4) cluster. Accordingly, we find 
three different phases depending on Ja^/J\\, as shown 
in Fig. [S] The term "phase transition" describes a 
recombination of VBS bonds. In the large-coupling 
regime J±/J\\ > 0(10), we can easily imagine that the 
on-rung pair is more stable than the on-leg pair and 
as many pairs as possible prefer to be formed on the 
rung bond [Fig. [2jV)]. The spin gap is therefore scaled 
with the binding energy of on-leg pair, i.e., A cx J||, 
which is consistent with the saturating behavior of A 
for J± ^ J|| . On the other hand, in the weak-coupling 
regime J±/J\\ < 0{1), all spin-singlet pairs are formed 
on the leg bond [Fig. mill)] because the binding 
energy of the on-leg pair is much larger than that of 
the on-rung pair. And, in the intermediate regime 
0(1) < J±/J\\ < 0(10), the spin-singlet pairs seem to 



be distributed in a balanced manner [Fig. [21JIV)]. Here, 
we notice an interesting relation to the phase transition 
in the S ~ ^ tube.— If we ignore a spin-singlet pair on 
each rung in the phases (IV) and (V) of our system, the 
remaining degrees of freedom are equivalent to those of 
the S = ^ tube. As it turns out, the recombination of 
the remaining VBS bonds between (IV) and (V) can 
be essentially equivalent to the phase transition in the 
5 = ^ tube. Then the (ignored) extra degrees of freedom 
yields the additional phase (III) in our S = ^ system. 

In conclusion, we study the ground-state properties of 
the three-leg S' = 3/2 Heisenberg tube with the DMRG 
method. It is confirmed that a spontaneous dimerization 
occurs and the spin-excitation gap opens for the whole 
coupling region. This may be a common feature of odd- 
leg half-integer-spin Heisenberg tube systems. We find 
that the gap for S* = | increases very slowly with increas- 
ing J± and it remains very small compared with J|| in the 
weak- to intermediate-coupling regions. For CsCrF4, the 
gap is estimated to be only a few K or less at normal 
pressures and, for example, additional condition such as 
applying pressure might be required to enlarge the ratio 
J±/J\\ in order to detect the gapped state. Moreover, by 
calculating the quantized Berry phase, it is shown that 
two phase transitions as recombination of VBS bonds oc- 
cur with varying the ratio J_l/ Jy although further work is 
desirable for quantitative evaluation of the critical points 
of the phase transitions. 
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